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Abstract 

The Hawking minisuperspace model (closed FRW geometry with a homo- 
geneous massive scalar field) provides a fairly non-trivial testing ground for 
fundamental problems in quantum cosmology. We provide evidence that the 
Wheeler-DeWitt equation admits a basis of solutions that is distinguished by 
analyticity properities in a large scale factor expansion. As a consequence, the 
space of solutions decomposes in a preferred way into two Hilbert spaces with 
positive and negative definite scalar product, respectively. These results may 
be viewed as a hint for a deeper significance of analyticity. If a similar struc- 
ture exists in full (non-minisuperspace) models as well, severe implications on 
the foundations of quantum cosmology are to be expected. 

Semiclassically, the elements of the preferred basis describe contracting 
and expanding universes with a prescribed value of the matter (scalar field) 
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energy. Half of the basis elements have previously been constructed by Hawk- 
ing and Page in a wormhole context, and they appear in a new light here. 
The technical tools to arrive at these conclusions are transformations of the 
wave function into several alternative representations that are based on the 
harmonic oscillator form of the matter energy operator, and that are called 
oscillator, energy and Fock representation. The framework defined by these 
may be of some help in analyzing the Wheeler-De Witt equation for other 
purposes as well. 
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1 Introduction 



The aim of this article is twofold: In the weaker sense it provides a formulation of 
minisuperspace quantum cosmology with a massive scalar field (the so-called Hawk- 
ing model) in terms of representations that are based on eigenstates of the matter 
energy operator E. Since this object does not commute with the operator defining 
the Wheeler-DeWitt equation, its eigenstates do not satisfy the latter, but their 
harmonic oscillator form motivates a change of representation for the wave function 
and an according transformation of the form of the Wheeler-DeWitt equation. All 
issues concerning the excitation of the scalar field oscillator modes become particu- 
larly transparent in the "energy representation" and related versions thereof. Also, 
our formulation may facilitate attempts to solve the exact Wheeler-DeWitt equation 
for whatever purpose. In the post-inflationary regime the matter energy eigenstates 
can be extended to approximate solutions of the Wheeler-DeWitt equation at the 
level of the WKB-approximation. 

In the stronger sense we present perspectives for the construction of exact wave 
functions that coincide with the approximate ones in the WKB-domain. Moreover, 
these states seem to have a preferred status, defined by analyticity properties in an 
expansion for large values of the scale factor. Although some of our conclusions are 
only conjectures, we seem to be able to define two exact Hilbert spaces of wave func- 
tions, playing a distinguished role. This is in some formal analogy to the one-particle 
Hilbert spaces of negative/positive frequencies in the flat Klein-Gordon equation. In 
this case the crucial property enabling us to decompose the space of solutions into 
two Hilbert spaces is Lorentz invariance. In contrast, the symmetries and covari- 
ance properties of the Wheeler-DeWitt equation do not suffice to provide such a 
decomposition fl|]. What we found is that the asymptotic analyticity structure of 
solutions may play a role analogous to Lorentz invariance. Accepting a preferred 
decomposition has of course implications for the conceptual basis of quantum cos- 
mology. What we cannot answer at the moment is the question to what extent this 
structure will apply for the full (non-minisuperspace) theory. 

The basic variables in the Hawking minisuperspace model 0] are the scale factor a 
of a closed Friedmann- Robert son- Walker universe and the value (f) of a homogeneous, 
minimally coupled massive scalar field. The external parameters are the mass m 
of the scalar field and a numerical constant p representing the operator ordering 
ambiguity in the Wheeler-DeWitt equation. Apart from some general remarks, we 
will work entirely within this model. 
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In Section 2, we introduce as tools for the analysis a representation in terms of the 
harmonic oscillator eigenf unctions associated with the matter energy, and a suitable 
Fock space notation. We have designed this Section as self-contained as possible and 
provided various formulae that are helpful in dealing with different representations 
of the wave function. For later reference, we distinguish between position, oscillator, 
energy and Fock representation. In Section 3, we use the fact that matter energy 
is approximately conserved after inflation to write down energy eigenstates that 
approximately satisfy the Wheeler-DeWitt equation in the corresponding domain 
of minisuperspace. These states have been used by other authors as well, at the 
level of the WKB-approximation. The asymptotic structure for values of a larger 
than classically allowed is related to the classical domain by an appropriate WKB 
matching procedure. The same asymptotic structure appears when the semiclassical 
WKB-expansion method is applied straightforwardly. 

In Section 4, we write down a fairly general expression for wave functions in the 
representation based on the original variables a and <fi. We impose an analyticity re- 
quirement in terms of an expansion in inverse powers of a that seems to single out an 
exact version of the approximate states considered in Section 3. In Section 5, a sim- 
ilar procedure, when applied to the Wheeler-DeWitt equation in the oscillator and 
Fock representations, seems to generate identical results. Assuming these to hold, 
we end up with a basis of exact and distinguished wave functions 5^ (a, <fi) , each 
describing (at the level of WKB-identification) an ensemble of collapsing/ expanding 
universes with matter energy (n + \)m. Half of these wave functions have been 
constructed previously in a wormhole context by Hawking and Page |§, and the 
way they emerge in our framework sheds new light on them. Using the indefinite 
Klein-Gordon type scalar product Q on the space of solutions of the Wheeler-DeWitt 
equation in Section 6, we decompose this space into two Hilbert spaces, with pos- 
itive/negative definite scalar product, respectively. If this decomposition is consid- 
ered as a preferred one, relevant for interpretation, positive probabilities may be 
written down by means of conventional Hilbert space techniques. As an example, 
we expand the no-boundary wave function in terms of this basis in Section 7. Such 
an expansion may be called "energy representation" in the sense of representing a 
state entirely in terms of expansion coefficients with respect to a basis of solutions 
of the Wheeler-DeWitt equation that describe universes of definite energy. Thereby, 
also the role of the above-mentioned Hawking-Page solutions as providing only half 
of a basis is illustrated. Some concluding remarks, concerning the tunnelling wave 
function and the significance of the structure we found for the conceptual issues of 
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quantum cosmology are given in Section 8. 

To conclude this introduction, we comment on the units used (see Ref. 0). Let 
a 2 = 2G/3ir. In what follows a tilde shall denote "true" physical quantities in units 
in which c = 1. The Planck mass and length are mp = (h/G) 1 ^ 2 and £p — h/rhp. 
The FRW space-time metric is given by 

ds 2 = -N{tfdt 2 + d{t) 2 dal (1.1) 

where da\ is the metric on the round unit three-sphere. The scale factor and lapse are 
rescaled as a = aa and N = aN. Let furthermore be the (spatially homogeneous) 
scalar field (with dimension h 1 ^ 2 ^ 1 , i denoting length) and m its mass, and set 
= cr7r\/20 and m = am. (A general scalar field potential would be redefined as 
V(<fi) = 2rr 2 a 4: V((j)). For the massive case we have V(i/>) = m 2 (j) 2 /2h 2 and V{4>) = 
m 2 cj) 2 / 2h 2 ) . According to this scheme we rescale the Planck mass and length as 
mp = amp = (2/i/37r) 1//2 and £p = a~ x tp = (3nh/2) 1 / 2 . The Lagrangian resulting 
from this ansatz, when expressed in terms of the rescaled variables, is thus 

r 1 ( ad 2 nT \ a 3 (ft Nm 2 l2 \ . . 

L= + Na)+—[- o-0 , (1.2) 

the action being S = J dtL, with dimension H. From now on, it is easy to restore the 
original variables at any stage of the quantization procedure. (This will be helpful 
for the reader who likes to go into the details of the semiclassical expansion carried 
out at the end of Section 3). By using units in which H = 1, the variable a as well as 
the mass parameter m become dimensionless, and we have mp = (2/37T) 1 / 2 w 0.46 
and £p = (3-7T/2) 1 / 2 w 2.17. The ratio between scalar field mass and Planck mass 
is thus m/mp = m/mp w 2.17m. In order to account for the necessary amount 
of density fluctuations fHH, we expect the mass parameter to be m ~ 10~ 6 , hence 
much smaller than 1. 

2 Representations for states 

Let us as preliminaries write down the equations governing the classical dynamics 
of the Hawking model MM, i.e. the trajectories (a(t) , <f>(t)) in the minisuperspace 
manifold {(a, 0)|a > 0}. If the space-time lapse N has been fixed (e.g. by assuming 
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a functional dependence N = N(a,(j))), the relation between the momenta and the 
time-derivatives of the minisuperspace variables is given by 



a da a 3 d(J) 

Pa = ~NTt P +=Ndt' (2 - 1} 

In the time gauge N — 1, the classical constraint equation reads 

a 2 + 1 = a 2 (0 2 + m 2 2 ). (2.2) 

When evaluated at some initial time to, it may be interpreted as a restriction on the 
set of initial conditions (a(t ), <j>(to), a(t ), <j>(to)). Once it is imposed for all times, 
the time evolution equations reduce to 

+ 3-0 + m 2 = O. (2.3) 
a 

The corresponding a-equation which arises from the Lagrangian formalism is auto- 
matically satisfied for all times on account of (|2.2| ). In a general time gauge, the 
above equations apply after replacing d/dt by iV -1 d/dt. 

We will not go into the details of the properties of classical trajectories. They 
have been studied by a number of authors (see e.g. Refs. @ [0). Let us just note that 
a typical (outgoing) classical trajectory leaves the inflationary domain ma|0| ^> 1, 
|0| 3> 1 when |0| settles to the order of unity. The subsequent evolution is of the 
matter dominated type, with undergoing rapid oscillations. Eventually (when the 
amplitude of falls into the domain of negative potential, ma |0| < 1), the scale 
factor a reaches its maximum value, and the universe recollapses again. During this 
process the energy of the scalar field 

E = y (0 2 + m 2 2 ) (2.4) 
is approximately conserved. The amplitude of the oscillations of is given by 

0ampl ~ m ■ (2.5) 

The maximum scale factor is a max w 2E. (In order not to deal too much with 
approximate quantities, one may simply define E = a max /2). The value of a at 
which the trajectory enters (and leaves) the domain |0| <1 (i.e. the value of the 
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scale factor at the end of inflation) is roughly given by a m i n ~ (E/m 2 ) 1 ^ 3 . Hence, 
only if a m i n <C a max , i.e. E ^> m -1 , we have a post-inflationary classical evolution 
at all. 

Canonical quantization is achieved by rewriting the constraint in Hamiltonian 
form and performing the usual substitutions (see e.g. ||). The result is the minisu- 
perspace Wheeler-DeWitt equation 

Hijj = 0, (2.6) 

where the state is represented as a wave function ip(a,<f)). Since in what follows 
different representations will be used, we distinguish between the notation of an 
"operator" and its "representation", in particular if derivatives with respect to a are 
involved. Let T> a be the operator that acts on a wave function ip(a, 4>) as the partial 
derivative d a . Then the Wheeler-DeWitt operator is given by 

H = V a V a + lv a + 2aE-a\ (2.7) 

a 

where 



represents the total matter energy (cf. £]4| for its classical analogue), and p is a pa- 
rameter accounting for the operator ordering ambiguity. There are good arguments 
in favour of p = 1 (see e.g. Ref. |§), but we will leave it unspecified. The last term 
—a 2 in ( p.7| ) represents the spatial curvature. In case of a spatially flat FRW model, 
one would omit it, and in case of an open FRW model one would change its sign. In 
analogy with ordinary quantum mechanics, the wave function tp(a,<f)) can be said 
to be in the position representation (which is just the defining representation here). 
To be explicit, the Wheeler-DeWitt equation in this representation reads 

(?#)(«, 0) = (d aa + -d a - \d H + m 2 a 4 (P 2 - a 2 ) ^(a, 0) = , (2.9) 

which is the standard form in which it is usually written down [0 and which deter- 
mines its form in all other representations. (We will not make attempts to modify it, 
e.g. by introducing different operator orderings or taking square roots of operators 
as inspired by the hope of making expressions simple in some particular representa- 
tion. In other words, we are not searching for an alternative wave equation, but just 
stick to (|2.9| ) as the starting point, although written down in other representations). 
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Let us begin our analysis by exploiting the fact that for any value of a the 
operator E (as acting on functions of 0) represents a quantum mechanical harmonic 
oscillator with frequency m and mass a 3 . Its eigenvalues are thus E n = (n + ~)m 
for non-negative integer n. Using the combination 

£ = m 1/2 a 3/2 0, (2.10) 

we define an alternative representation of states by 

ip(a,<j>) = m 1/4 a 3/4 ^(a,0 . (2.11) 

Since £ plays the role of an oscillator variable, one could call ^(a, £) to be in the 
oscillator representation. The matter energy operator becomes essentially a unit 
harmonic oscillator, 

£ = -(-%+a. (2.12) 

The operator T> a (which was d a in the representation ijj(a, 4>)) takes a different form 
now. Let E) a be the operator that acts on a wave function in the representation 
■0(a,£) as the partial derivative d a . Then we have 

V a = E) a + ^-IC (2.13) 

with 

K = \ {0,94 = 0^ + I = = + I. (2.14) 

The Wheeler-DeWitt operator in the oscillator representation ip(a,£) is now still 
given by fl2.7j) , but with D a being represented as fl2.13|) , and E) a being represented 
as d a . Thus, one may write 

H = ID a lD a + P -lD a + ^-lD a + K + ^ /C 2 + 2a£ - a 2 (2. 15) 

which is valid in any representation (just as (|2.7| ) is), as long as by -D a the appro- 
priate operator representation is understood. 

The eigenfunctions of E are just those of the unit harmonic oscillator with co- 
ordinate £. In terms of Hermite polynomials (which are generated by e 2 ^*~* = 
En=o ^ n H n (t)) they read 
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with n a non-negative integer. By expansion with respect to these, we define a 
further way of writing wave functions 



feO = £/*( a )*»(0, (2.17) 

71=0 

where the component functions f n (a) may be regarded as providing the state in 
the energy representation (although we will encounter a further meaning of this 
word later on). This notation is justified by the fact that the operator E is now 
diagonal: it sends f n {a) to E n f n (a). Its action may symbolically be written as 
(Ef) n (a) = E n f n (a). Note that even n belongs to the even (if) (a, — (f>) = ^( a )0)) 
and odd n to the odd (if} (a, —(f)) = —if>(a, 0)) sector of wave functions. Since the \&„ 
are an orthonormal basis, ( [2.17 ) may be inverted to give the oscillator excitations 



/oo 
df* n (0V(a,0 (2-18) 
-oo 



in terms of the wave function in the oscillator representation. 

In performing ( |2.17| ) we have implicitly assumed that the wave function if) is 
sufficiently well-behaved for large £ (or <f>) so as to allow for such an expansion. A 
quite restrictive condition on general wave functions would be square integrability in 
the matter variable £ (or 0, which is equivalent). Although this would offer a Hilbert 
space structure for any value of a, the more interesting candidate wave functions of 
the universe are of distributional character with respect to this structure, and we 
just assume that the expansion (|2.17|) is possible. The formal squared norm of wave 



functions 

d<f>r(aA)^A) = / dZPM^M = £/>)/» (2-19) 



n=0 



may but need not be finite. The prefactor m 1 / 4 a 3//4 in (|2.11|) has been chosen so 
as to give these expressions a simple form. The formal hermiticity of operators like 
id<p, id^, ilC and E with respect to the according scalar product is evident. 

The Wheeler-DeWitt operator in this representation is given by ( [2.15D with E) a = 
d a , E acting diagonal with eigenvalues E n and JC acting as (in a symbolic notation) 

(Kf) n (a) = - y n (n-l)f n -. 2 (a) + \j(n + l)(n + 2) f n+2 (a) , (2.20) 
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its square being given by 

(/C 2 /)n(a) = y n (n-l)(n-2)(n-3)f n ^(a)- 

l(n( n - 1) + (n + l)(n + 2)) / n (a) + (2.21) 

^(ra + l)(n + 2)(n + 3)(n + 4) / n+4 (a) . 

For completeness, we write down the Wheeler-DeWitt equation in this representation 
explicitly: 

(Hf) n (a) = (d aa + V - d a )f n (a) + - a d a {Kf) n {a) + l) (/C/) n (a) + 

^ (/C 2 /)n(a) + (2a£ n - a 2 ) / n (a) = . (2.22) 

It involves differences with respect to n rather than derivatives. Note that the 
highest component function f n +±{a) appears algebraically. Hence, the Wheeler- 
DeWitt equation in the energy representation is just a recursive expression for f n (a) 
(n > 4) in terms of (fo(a), /3(a)). 

The form ( | . 1 2|) of E amounts to define formal annihilation and creation opera- 
tors 

* = ^fi*-#zm a * = &- Bd (2 - 24) 



from which 



and 



2 

The square of the latter turns out to be 



E = m (j$A +l)= m { Af+ \) ( 2 - 25 ) 
K = I (A 2 - (A^) 2 ) . (2.26) 



j(^ 4 -AT(Ar-l)-(AT+l)(Ar + 2) + (^ t ) 4 ) . (2.27) 
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The index n of the components / n (a) represents the eigenvalues of the oscillator 
number operator M = A^A. The dagger denotes hermitean conjugation in a formal 
sense, with respect to the scalar product associated with the squared norm expres- 
sions (|2.19|) . By formally identifying ^/ n (0 with the abstract state |n) and using the 
commutator relation [A, A^] = 1, we find the usual structure defining a Fock space 



A\n) = y/n \n — 1) A ] \n) = Vra + l|ra + l>. (2.28) 

The n-the eigenstate is generated out of the ground state (note that *A|0) = 0) by 

|n) = |0) . (2.29) 

vnl 



A given wave function may be written in the form 

KM = J2f n (a)\n) = ^f n {a)^\0)=F(a,Ai)\0), (2.30) 

n=0 n=0 V n - 



thus defining a "Fock representation" for states. If (if), a\ is defined as J2n{ n \fn( a )i or 
equivalently as (0|.F*(a, A), the formal squared norm (|2.19|) is given by (if),a\i(),a). 
The analogue of ( |2.18| ) is 

f n (a) = (n\^,a), (2.31) 

and the orthonormality of the oscillator basis carries over to (r\s) = 5 rs . If \ip,a) 
is written as J 7 (a, A^) |0), the operator A is formally represented as d/dA^ , whereas 
A^ may be considered as a multiplication operator, and ID a is the partial derivative 
d a . The relation between ip(a, £) and jF(a, A') has so far been given only through a 
number of intermediate steps, involving Hermite polynomials and an infinite sum. 
It can be made more explicit by noting the identifications of functions of £ with 
abstract Fock space states 

oo 

e^ = ]T« n (A;)^ n (0 < — ► tffy/2e-$*j V * kA * + 3W\0), (2.32) 

n=0 

for fixed k and 

OO 1 

5(t-ri) = £tt„fa)*»(0 ^^e-^e^'-K^lo) (2.33) 



n=0 - /7r 



for fixed 77. A further interesting relation illustrating the appearence of Hermite 
polynomials is 

C< — - <f^V2(-i) n 2~ n l 2 H n (iA ] )e^ Ar)2 \$) . (2.34) 



10 



for non-negative integer n. Using ( |2.33| ), one finds for a given wave function if}(a,£) 
that 

\/7T J-oo 



Inverting this relation is not that straightforward. Defining if} (a, k) to be the inverse 
Fourier transform of if)(a,£), i.e. 

^ POO rjh _ 

iP(a,0= ^=iP(a,k)e^ } (2.36) 

J—oo v Zix 



and making use of (|2.32j ), equation ( |2.35| ) may be expressed alternatively as 

dk 



T{a, J$) = <f^V2 [°° if>{a, k)e~ h^-i^f & - . ( 2 .37) 

J — oo \/2,TT 



This may be inverted to give 



$ a> A,) = _L r ^t )e - i^t+ws*)^- 4 *» (2 38) 

■^7T J-oo V27T 



Combining this last equation with ( |2.36D and interchanging the A* and k integrations 
gives a direct formula for ifi(a,£) in terms of J-"(a,A'), but with an integrand that 
exists only in a distributional sense. 

The Wheeler-DeWitt operator in this representation is given by Q2.15| ) with E) a = 
d a , E from ( ^25| ) and K from (ggg) . Note that, due to (|2T27|) , the operator /C 2 
contains ^4 4 . This makes the Wheeler-DeWitt equation in the Fock representation 
a fourth order differential equation for !F(a,A'). 

We are thus able to represent states in various forms, namely as if} (a, <f>), if} (a, £), 
f n (a), \if>, a) and jF(a, A^). The latter three forms are of course closely related to each 
other, referring to the scalar field energy as a variable. The different versions of the 
Wheeler-DeWitt equation are equivalent, which makes the choice of representation 
a matter of convenience and taste. 



3 Approximate solutions 

The framework of representations as given in the last Section is "kinematic" in 
nature, i.e. it makes no reference to the Wheeler-DeWitt equation. It is clear that 
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an object like |n), for some fixed n, does not even approximately satisfy it. On 
the other hand, we know that classically the energy E is approximately conserved 
in the post-inflationary regime. We thus expect that there are solutions of the 
Wheeler-DeWitt equation which behave to some accuracy like 

hM) = /i(a)|n) (3.1) 

for values of a which belong to the post-inflationary classical domain (i.e. between 
end of inflation at a m i n »s (E^/m 2 ) 1 / 3 and maximum size at a max ~ 2E n ). 



Inserting (|3.1|) into the Wheeler-DeWitt equation in the Fock representation of 
( 2.15 ), the reason why it cannot be an exact solution turns out to be that K, and K? 



as represented by (|2.26|) and (|2.27|) contain powers of A and A^ that mix the elements 



|r) of the oscillator basis. These terms arise from the a-dependence of E in 
When a is large, they can be neglected, and the harmonic oscillator dynamics follows 
the dynamics of the gravitational sector. This leads to the adiabatic approximation, 
a technique which is frequently applied (see e.g. Ref. [0]). In order to have a non- 
trivial classical domain at all (a m i n -C a max ) we must choose n ^> vrC 2 « 10 12 , hence 
much larger than 1. Neclecting all non-trivial powers of A and A* in (|2.15 ), we pick 
up the effective potential term — ^{n 2 + n + 1) from JC 2 . Taking into account the 
operator ordering term as well, we end up with an approximate equation for the 
range a > a min 

(d IM + -d a + U(a))h(a) = (3.2) 

with 

U(a) = 2aE n -a 2 - -^(n 2 + n + l) . (3.3) 
8a 2 

From the point of view of the oscillator basis, this is just the projection of the 
Wheeler-DeWitt equation with (|3.1|) inserted onto the n-the basis element. Hence, 
(|3.2|) is identical to (n\Hh(a)\n) = 0. It is the best we can do when all oscillators 



other than \n) are ignored. 

By redefining h(a) = a~ p ^ 2 g(a), one kills the first order derivative. In the range 
a ^> a m j n , the a -2 contribution to U(a) (together with an additional term from the 
above redefinition) is small. Neglecting it, the new approximate equation reads 

(d aa + 2aE n -a 2 )g{a) = 0. (3.4) 

The solutions thereof may be expressed in terms of parabolic cylinder functions 
[|T0|j. Since this is not very instructive, let us look at some possible ranges of a. For 
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a <C a max , the potential in ( |3.4J) is dominated by 2aE n , the according solutions being 
the two Airy functions Ai(— [2E n ) 1 ^ a) and Bi(— (2£' ri ) 1//3 a). Since the arguments 
— (2E n ) 1 ' 3 a of these functions are much less than —1, we are in the range in which 
they oscillate rapidly. Up to a multiplicative constant, the according asymptotic 
expressions for g(a) are 

cT^oscfJ (2£ n )V2 a 3/2 + £) (3 5) 

with osc = cos for Bi and osc = sin for Ai. 

When a approaches a max , the —a 2 contribution in the potential becomes im- 
portant and will slow down the oscillations. Here we note that so far it was not 
necessary to restrict a to be less than a max . At a ~ a max « 2E n , the potential 
vanishes. For a ^> a max , the overall behaviour is exponential, and we find two 
asymptotic solutions for g(a) of the type 

-l/2=Fi£/2 exp f ± 1 2 T ^ (3 6) 

For given £7 n , this range is "classically forbiden", because all classical universes with 
matter energy E n cannot extend to such sizes. Note that for n<m~ 2 , the asymptotic 
solutions (|3.6|) still exist, athough they do not have an oscillatory domain of the type 



So far we have not stated anything about the allowed range of <fi. The classical 
amplitude of oscillations is given by Q2.5|) , and we expect the approximation for the 
oscillating regime ( |3.5| ) to hold if |0|^0ampi- In terms of £, this means \^\<\/2n. 
This is in fact just the domain in which the oscillator basis ^ n (0 is considerably 
non-zero. (Note that the amplitude of the classical oscillations with energy mn 
corresponding to the Hamiltonian (|2.12| ) is just £ amp i = v2n). In the exponential 
regime (|3.6|), an identical condition on £ holds. Hence, for a ^> a m \ n the interesting 
values are concentrated in a stripe around the = axis that gets arbitrarily narrow 
as a increases. 



In order to link the asymptotic form ( |3.6|) with the oscillating behaviour ( |3.5| ), 
we may apply the standard WKB matching procedure between domains in which 



the potential has different sign [Tl]. Introducing the number 



©»= + ( 3 - 7 ) 
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we fix the normalization in the oscillatory region by defining two real approximate 
solutions of the type ( |3.5|) 



a) = a-^- 1 " (2E n )- 1/4 cos (e„ - ~ {2E n ) 1 ' 2 a^ \n) (3.8) 
\Q-,a) = a~^- l l' (2E n )-^ sin ( 9 n - ~ {2E n )^a^\ \n) . (3.9) 



The result of applying the WKB matching procedure (which is too boring to be 
shown here) is that, in the exponential domain, 

/ On \ ~ E n/ 2 / 1 1 \ 

|n+a) = a- p/2 - 1/2 (^) exp (±o a - E n a + ^) \n) (3.10) 



|n~ a) = ^" P/2 - 1/2 (^) n exp(-ia 2 + E n a-l^)|n). (3.11) 

For n<m~ 2 there is no oscillatory region, and we define the asymptotic normalization 
by the these two expressions as well. The factor \ in (|3.11|) may look a bit stange, 



but it is an immediate consequence of the matching procedure (it is actually part of 
the standard formulae) and the fact that we have arranged the normalization of the 
oscillating behaviour ( |3.8| )-( PT9"| ) in a symmetric way. Note that the Airy functions 
Ai(x) and Bi(x), when expanded for x — > — oo and x — > oo (i.e. in the oscillatory 
and in the exponential domain) display a precisely analogous factor | for Ai(rr) (cf. 
also (|3J)). 

Defining two other sets of approximate solutions 5^ by 



or, inversely, 



S± = e ±ie "(n+ T zn;), (3.12) 



K = ^(e-* "2+ + e^~;) (3.13) 
K = \ (e- 40 " - e i0 « 2") , (3.14) 



we find, in the oscillatory regime, 



|S±> = a-^ l '\2E n )~ l " exp (±| (2E n )^a^ ) \n) . (3.15) 
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The involved structure of prefactors was necessary in order to achieve this simple 
form. The significance of the normalization of this expression (in particular the 
factor {2E n )~ l l A ) will become clear later. 

Thus, we have a set of approximate solutions at hand that correspond to pre- 
scribed values of the matter energy. These functions have been used by Kiefer in 
his discussion of wave packets in the Hawking model 0. Also, the exponential 
behaviour of f2o (a, <f>) has been displayed by Page || as the dominant part of the 
no-boundary wave function for small and large a. The Q~ also occur in the work 
of Hawking and Page |§ , to which we will refer later on. 

The physical significance of any of the wave functions with n 3> m~ 2 at the 
WKB-level is to represent an ensemble of contracting/expanding classical universes 
with post-inflationary matter energy E n . The expressions ( |3.15| ) may be viewed as 
semiclassical WKB-states around a family of expanding classical backgrounds with 
action S (a) = - \\f2Ea il2 (th is is discussed in a bit more detail in Ref. fOl). 
Note however that the classical matter energy E is approximately conserved only 
in the domain \(f)\ <C 1, i.e. "after" inflation. Hence, the individual |S^,a) can 
be expected to contain essentially a single oscillator excitation (i.e. of \n)) only 
if a ^> a min « (En/m 2 ) 1 ^ 3 . For smaller values of a we expect |H*, o) to be a 
non-trivial superposition of (virtually) all oscillators |r) and thus all approximate 
expressions we gave for these wave functions to break down. Furthermore it is likely 
that these states are not well-behaved for small a. This is because there are many 
classical trajectories to some given value of a max which behave quite singular as 
a — > (namely of the "collapse" type \<f>\ — > oo). Also, we cannot expect |S^,a) 
to satisfy nice properties near the zero potential curve a 2 m 2 (p 2 = 1, along which 
usually "nucleation" is assumed to occur. In the case n<m~ 2 we can talk about 
pure tunnelling states that do not correspond to classical universes at all. 

Although, so far, the wave functions have only been identified as approximate 
solutions of the Wheeler-DeWitt equation, it is reasonable to suppose the existence 
of exact solutions that behave qualitatively in the same way. In the range a ^ ct m i n 
this means that the excitations of oscillators other than \n) may be non-zero but are 
small. Moreover, we expect to obtain a basis for the set of solutions: At some fixed 
value of a, the degrees of freedom contained in {fi^} correspond to two free functions. 
Hence, at the approximate level, any initial data (ip, d a ip)\ a=ainl may be expanded in 
terms of the oscillators \n) (i.e. in terms of the functions \I/ ra (m 1//2 a 3 / 2 0) = \I/ n (£)). 
This behaviour is expected to carry over to the exact case. 



15 



For large a, the oscillator mixing operators (powers of A and A') contained in 
K, and K? get suppressed. Hence we expect the approximation of the predominance 
of one oscillator \n) in a wave function to become arbitrarily accurate as a — > oo. In 
this limit, we may treat all states belonging to non-negative integers n at the same 
footing (including the pure tunnelling states for small n). One must of course be 
aware that an exactification of our wave functions will leave a considerable amount 
of freedom. However, by combining the large-a expansion with the assumption of 
certain analyticity properties, we will encounter an appearently distinguished way 
to single out a unique exact solution for any choice of the ± label and any n. 

Let us at the end of this Section provide another argument leading to the same 
exponential behaviour of wave functions. One may apply the standard semiclassi- 
cal WKB-expansion scheme [[TB[], based on a Born-Oppenheimer approximation, by 
rewriting the Wheeler-DeWitt equation ( |2.9| ) in terms of units that make the Planck 
length explicit. At the end of the introduction we have displayed the relation be- 
tween the "true" quantities a, 0, m and the rescaled ones a, and m. Restoring 
the true units, one finds the kinetic part of the gravitational field multiplied by 
the Planck length squared and the curvature contribution multiplied by the inverse 
of the Planck length squared. The matter sector remains unaffected (up to a nu- 
merical factor of order unity). This structure may equivalently be written down 
by introducing a formal book-keeping parameter A (playing the role of the Planck 
length) that is treated as a small quantity (and reset equal to 1 in the end). The 
Wheeler-DeWitt equation (|2.9| ) thus becomes 

\ 2 (d aa + -da)-\ fy* + ™V</> 2 - ^ ] V(o, 0) = . (3.16) 

Ob a /\ j 

This may formally be achieved by replacing a — > X~ x a, — > A0 and m — > Am in 
( |2.9| ). According to the WKB-description we expand a wave function as 

^(a, 0) = exp (^{^S. + Sl (a, 0) + X 2 S 2 (a } 0) + A 4 S 3 (a, 0) + . . . )) . (3.17) 

This is the direct ("naive") semiclassical treatment of the Hawking model. It essen- 
tially assumes that the gravitational field variable a is " (quasi) classical" , whereas 
the matter field plays the full quantum role. Inserting the ansatz (|3.17|) into (|3.16| ) 
and isolating powers of A 2 yields a sequence of equations for the Sj(a,(f)). In the 
end, we restore the original Wheeler-DeWitt equation by setting A = 1. 
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The equation at 0(A~ 4 ) turns out to be just c^So(a, <f>) = 0, telling us that 
So = Sq(o). At 0(A -2 ) we obtain the equation (d a So(a)) 2 = —a 2 . Formally, 
this is the Hamilton- Jacobi equation for the pure gravitational field and thus fits 
into the general semiclassical scheme. On the other hand, it has no real solu- 
tion, which just reflects the fact that an empty closed FRW universe does not 
exist. Proceeding straightforwardly, we write down the two imaginary solutions 
So{a) = ^ia 2 /2. The WKB-phase factor e lSo ^ = e ±c? l 2 thus provides already the 
dominant large-a behaviour of Q^(a,(p). Rescaling the next order contribution as 
e iSi{a,4>) _ a -i/2-p/2 x ( aj ^ we find at o(X°) the Wheeler-DeWitt equation to state 

(±d a + E) X (a,<l>) = (3.18) 
with E from fl2.8|) . This is the step where usually the effective Schrodinger equa- 



tion (as a minisuperspace version of the Tomonaga-Schwinger equation) arises |13 
Obviously, we encounter a Euclidean type Schrodinger equation. In the adiabatic 
approximation, one neglects the a-dependence of E. By introducing the oscillator 



variable £ from ( |2.10| ), the energy operator becomes represented as ( 2.12 ). The adi- 



abatic approximation amounts to keep the derivative d a unchanged in ( |3.18|) . Thus, 
factorizing x = m _1//4 a _3//4 x (cf. fl2.11|) ) into a product A(a)X(£) (in the usual 
context this would be a stationary state) yields immediately an eigenvalue equation, 
hence X(£) — *& n (0 f° r some non- negative integer n, as well as the "time evolution" 
prefactor A(a) = e TEnCl . Putting everything together, we have reproduced qualita- 
tively the behaviour ( QD )-( STTT ) of the wave functions fl^ . The additional factors 



of the type a TE ™/ 2 are effects at 0{\ 2 ). 

Thus, for large a, the approximate solutions (a, 0) may be considered as 
semiclassical states, built around a "pure tunneling" background gravitational field 
(which nevertheless — in the way how it appears in a WKB-expansion — displays 
formal similarities to a true classical background). This way of looking at things 
may seem a bit strange, but since we are faced with a wave equation (as opposed to 
classical equations of motion) we cannot exclude that the tunneling region a ^> a max 
plays an important role in the structure of the space of solutions or in the conceptual 
foundations of a quantum theory of the universe. This may be in some correspon- 



dence with the idea of a Euclidean path-integral [^4[||15|| or some other principle 
which provides some additional structure that is invisible for the semiclassical tech- 
niques. Anyway, it may be taken as a further motivation in favour of examining the 
limit a — > oo of the Wheeler-DeWitt equation. 
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4 Exact solutions in the position representation 



Leaving the level of approximate wave functions, we will try now to define a large-a 
expansion scheme that enables us to specify a set of exact solutions of the Wheeler- 
DeWitt equation. An appropriate ansatz is modelled according to the asymptotic 
behaviour of as given by (|3.6|) , multiplied by \n), i.e. by the function \& n (0- The 
dominant behavour for large a and fixed of the latter is given by the exponential 
e -C 2 /2 = e -ma 3 4> 2 /2 ^ gee e q Ua ti ons (|2.10|) and (|2.16|) ). This is followed by the factor 



e ±a /2_ fj ence; j n a large-a expansion at fixed 0, the leading order in the exponen- 
tial is a 3 , with a 0-dependent coefficient. Separating the positive powers of a, as 
appearing in the exponential, from a pure a q term (which may represent e 9lna ) and 
expanding the remainder in terms of a -1 , we write down as a general ansatz in the 
position representation 

V(a, 0) = a q exp (F 3 (0)a 3 + F 2 (0)a 2 + F^a) £ . (4.1) 

This looks fairly general, the only severe restriction being that only one exponential 
of the above type is involved. In a paper on wormhole solutions, Hawking and Page 
U have effectively used a similar type of expansion, and some of the structure we will 
encounter appears there as well. By inserting the ansatz into the Wheeler-DeWitt 
equation ( |2.9|) and separating powers of a, one obtains a sequence of ordinary differ- 
ential equations for the functions F r (<f>) and G>(0), thereby expecting the freedom 
of choosing arbitrary integration constants. 

It is important to note that in an expansion in negative powers of a, the in- 
dependent coefficient functions cannot in general be expected to display regular 
behaviour. This may be illustrated by expanding the function (1 + a 4 4 ) 1 ^ 2 (which 
is C°° on the domain a > 0, = arbitrary) as a 2 2 + |a _2 -2 + . . .. For fixed 
0, the domain of convergence of this series is a > |0| _1 , which breaks down at 
= 0. The reason for limited convergence is that the function has singularities at 
a = ±\/50 -1 , against which a condition like C°° is insensitive. 

However, we are free to impose analyticity properties on the coefficient functions 
in the ansatz for if} (a, 0), as long as they are compatible with the Wheeler-DeWitt 
equation. By doing so, we touch upon a deeper structure which is not yet completely 
understood. We choose as a first condition the most natural requirement: 

(i) The functions F r ((j)) andG r ((p) are real analytic, i.e. they are real and analytic 
in a neighbourhood of the real axis in the compex <p-plane. As a consequence, they 
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admit a Taylor expansion at = 0. 

Hence, by truncating the series in ( |4. 1|) at fixed and sufficiently large a one should 
obtain a good numerical approximation for ip. There will be a function i?(0) such 
that the series converges for all a > -R(0). This defines a domain of convergence in 
minisuperspace. In order to exclude catastrophic behaviour of the leading a 3 term 
in the exponent, we require in addition 

(ii) F 3 (0) < for all 0. 
The significance of this condition may be illustrated for the case of the approximate 
expressions for Vt^- It allows for the dominant term I 2 = e~ ma * I 2 which appears 
in ty n (f;), whereas it prevents a behaviour of the type I 2 = e ma * I 2 which would 
otherwise be possible. 



Remarkably, the ansatz (|4.1|) together with the two conditions (i) and (ii) seem 



to leave only a discrete freedom. We believe that the general solution is characterized 
by a non-negative integer n and a choice of sign, i.e. a two-valued label ±. The 
wave functions appearing in this way are reckognized as exactifications of f2^(a,0) 
(and henceforth called by the same name) that are distingushed by their analyticity 
properties. 

Let us look at the sequence of equations generated by inserting ( |4.1|) into 
and dividing by (|4. 1|) . At the leading order 0(a 4 ) we find 



F^<p) 2 = 9F 3 (<f ) ) 2 + m 2 <f ) 2 . (4.2) 

Condition (ii) specifies a unique solution. Let us suppose that -F 3 (0) ^ 0. As a 
consequence, F 3 (<fi), which is given by a square root whose argument never vanishes, 
will always be non-zero. For — > ±oo it will either tend to oo or — oo. In both 
cases it is not possible for F 3 (0) to be non-positive for all 0. Hence, we must have 
F 3 (0) = 0. This fixes F 3 (0) = 0. Further differentiation of (|4.2j ) generates the 
expansion 

F 3W = -m (I ^ + 1^+1^ +...), (4.3) 

condition (ii) fixing the sign. This also agrees to leading order with the behaviour 
of the exponential factor e - ^ I 2 = e~ ma ^ I 2 in \l/ n (£). The behaviour of F 3 ((p) for 
large is 

F 3 (0) w -0.047 me 3 *, (4.4) 
where the constant has been determined by numerical methods. 
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At the next order 0(a 3 ) we find the equation 

F^(0)^(0) = 6F 2 (0)F 3 (0) . (4.5) 

Since F 3 ((p) is already uniquely determined, -F 2 (0) is fixed up to a multiplicative 
constant. In particular, we find i^(0) = 0. The equation at 0(a 2 ), reads 

2i^(0)i^(0) + i^(0) 2 = 6^(0)^(0) + 4F 2 (0) 2 - 1 . (4.6) 

Again, -Fi(0) is determined only up to an integration constant. Inserting = 0, we 
get 4F 2 (0) 2 = 1. Thus, there are two possibilities ^(0) = ±~. For either sign, the 
solution of fl4.5| ) is now uniquely determined to be 

F 2 (0) = ±fi + -0 2 + — 4 + — 6 + ..) . (4.7) 
V2 4 r 128 r 512 r / v ' 

The leading behaviour of .F 2 (0) for small thus reproduces e ±a2/ ' 2 , which we already 
know from (|3.6|). For large we find 



F 2 (0) « ±0.231 e 2 ^ (4.8) 

as the leading asymptotic behaviour. 

At O(a), the number g as well as the operator ordering parameter p and the first 
function of the series in (|4.1|) come into play. The equation reads 

Fgffl + 2i^(0) H| + 2^(0)^(0) = 4F X (0)F 2 (0) + 3(2 + p + 2g)F 3 (0) . (4.9) 

Since, by condition (i), Gq(4>) is analytic, its Taylor expansion at = exists and 
starts with ki<p n + k 2 <p n+l + ■■■■, where n is some non-negative integer and k\ 7^ 0. In 
the limit — > we have F^(<f))G' ((j)) / Go(4>) — > —mn, and Q4.9|) reduces in this limit 
to -Fi(O) = =f(?t, + |)m = =F-E„. As a consequence, equation ( |4.6j ) admits a unique 
solution 

F x (0) = T E n + -J- ( - 1 T 4mE„)0 2 + -J- (3 =F 2m£ n )0 4 + 
lom v ' 25om v y 

5j l_(41 T 4 8n , B „)/ + ... (4.10) 

The leading term thus reproduces e =F " Ena from ( |3.6|) . For large we find 

Fi(0) « ^0.654 £ n e (4.11) 
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as the dominant behaviour. 

From now on we just describe the general structure of the subsequent steps. The 
limit <p — > of the 0(a°) equation yields that <f)~ n Gi(<j)) must be regular. Moreover, 
the leading order part of this equation fixes 

q = - + -- P -T-El (4.12) 
2 42 2™ v; 

The first term 3n/2 combines together with <p n from Go(0) an d a m™/ 2 from the 
overall normalization into £ n — cf. ( |2.10|) — , which provides the leading order of 



the n-th Hermite polynomial H n {£). Taking into account the factor a 3 / 4 between 
the position and oscillator representation (see ( |2.11| )), we exactly reproduce the 
contribution a~ p / 2_1 / 2=F E ^l 2 from (|3.10|) - (|3.11|) . This result is inserted into ( [4.9|) , by 



which Go(4>) becomes unique up to the multiplicative constant k\ which survives as 
an overall normalization freedom for r/>(a, <fi), and we find 

G (<f>) = k x r (l + ^- (3mE n -m 2 T lT 2m 2 E 2 n )<f ) 2 + ...). (4.13) 

The overall pattern seems to persist at all orders. All equations are of the linear 
inhomogeneous type, leaving an integration constant which is determined at the 
next order by the analyticity requirement. We have checked this up to Gf {<))). 

By re-writing the series in as an exponential, one may re-arrange terms 

in a more explicit way. This is in fact what Hawking and Page have done for 
the case p — 1. When translated to our formulation, their result seems to make 
explicit how the pattern determining the functions G r {<p) persists to all orders. The 
uniqueness of the coefficient functions is not considered as an important issue in 
Ref. ||, but it is effectively achieved by throwing away ln(0)-terms at each order). 
Since these authors intended to construct wormhole solutions, they considered only 
the exponentially decreasing behaviour e~ a I 2 . They arrive at a set of functions 
\l/ n (a, 0), which we will denote as \P^ p (a, 0). 

Hence, without having a rigorous proof, we conjecture that all G r (<p) exist and 
are uniquely determined, once the ± label and n have been chosen. The case of the 
exponentially decreasing sector (lower sign) for p = 1 seems to be covered by Ref. 
[H]. Moreover, recalling the physical discussion of the approximate wave functions 
in Section 3, we believe that the series Q4.1D defines exact solutions fi^(a,0) of the 
Wheeler-DeWitt equation for all (a, 0), i.e. also outside the domain of convergence 
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of the series (in case this domain does not agree with the whole of minisuperspace) . 
We also note that, pulling an overall factor n out of f2^(a, 0), only even powers of 
remain, hence f^(a, — 0) = (— ) n f2^(a, 0). The functions \I/^ p (a,0) as displayed 
by Hawking and Page appearently coincide (up to normalization) with our fi~(a, 0). 
(The leading order of their equation (72) is just the explicit formula for the n-th 
Hermite polynomial). 

Due to the expansions of the coefficient functions around = we have an 
idea how f^(a,0) behaves for small and large a. By looking at the expression 
(fO|) of -F 3 (0) for large 0, we expect that, for sufficiently large and fixed a, the 
dominant behaviour for large is exp(— 0.047 ma 3 e 3r ^). If this conclusion holds, the 
wave functions are actually more dampted than one would expect from the factor 

/ 2 = e~ ma ^ I 2 of the oscillator basis fl2.16p alone. One may in fact perform an 
independent analysis of solutions of the Wheeler-DeWitt equation in terms of large 
0, thus identifying the fl^(a, 0) by means of the asymptotic expressions ( [4.4] ), ( |4.8| ) 
and ( |4.11| ). We will not go into these details but just complete our conjecture by 



noting that our wave functions are likely to be well-behaved as |0| — > oo. 



5 Exact solutions in the oscillator and Fock rep- 
resentations 

So far we have considered an expansion in a -1 at constant 0. For actual computa- 
tions the oscillator and Fock representations turn out to be more suitable. As we 
have seen, the series in (|4.1j) contains an overall factor n which, together with a 3ri//2 
from ( [4.12] ) makes up an overall factor £ n . Since the remainder contains only even 



powers of 0, the transformation to the variable £ — cf. Q2.10 ) — amounts to substi- 
tute 2 — > m _1 a~ 3 £ 2 , and even powers thereof. As a consequence, we never pick up 
half-integer powers of a or negative powers of £. Moreover, keeping £ fixed means to 
follow a curve in minisuperspace whose 0-coordinate value decreases towards the axis 
= as a — > oo. This should not affect questions of existence and convergence very 
much (or even improve the situation). The finite sum F 3 (0)a 3 + F 2 (0)a 2 + Fi(<p)a 
reduces to the expression ±a 2 /2 ± E n a plus a series of functions that contains only 
negative integer powers of a. 

We can thus re- arrange the expression (|4.1| ) in terms of a and £. In order not to 
overcomplicate things, we assume the choice of the ± sector and of n has already 
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been made, and explicitly insert the leading orders. Transforming if) (a, <f>) into the 
oscillator representation if)(a,£) by (|2.11|) , we end up with 



nJ (a, = Cn ^ exp (± 1 -a 2 T E n a) £ ^1 , (5.1) 

where the functions are analytic and is an arbitrary overall normalization 

constant. The numbers q n may either be infered from Q4.12Q by taking into account 
the correct transformation factors, yielding 

£ = -\-\t\eI, (5.2) 

or left undetermined in order to be re-discovered in the oscillator or Fock formalism. 
In view of our conjecture, the should be uniquely determined. Note however 

that these functions arise from a rearrangement of orders in the series of ( f4.1| ), 
combined with a series stemming from the exponent. Computationally, they are 
related to F r (<f>) and G>(0) in a non-trivial way. Due to the discussion given in 
Section 3, we expect the eigenfunctions of E to appear at leading order in a -1 . 
Hence, we set 

i^(0 = Vni« r »(e), (5-3) 

which will be justified simply by being consistent. The prefactor \fn). is for later 
convenience. 

This formulation is still a bit awkward for general n. We just report briefly on 
the expansion for n — 0. Redefining the sum J2 r ^ro(£)/° r as ex P(Xr 9t{0/ ar )-> one 
finds that all <7^(£) are polynomials in £ containing only even powers and being of 
order |(r + 3 — j) with j — 0, 1 or 2. The formal criterion necessary to single out this 
unique sequence of functions turns out to be the exclusion of homogeneous solutions 
for the g^(£) of the error function type. This just prevents a behaviour in the wave 
function that would contradict condition (ii) of Section 4. 

Some simplification occurs by translating the above expression OOP into the 
Fock representation in which states are written as J-(a, A^)\0), and this is the setup 
we will consider now in more detail. The function Q is just (A*) n \0). All other 
F±(0 - - which are of the form Gf n {A ] )\ti) - - are written as (A^) n Gf n (A*)\0), 
thereby defining a set of functions Gf n (A^). The wave functions thus become 

a) = C± a«* exp (± l -a 2 T E n a) (A*) n £ , (5.4) 



=0 
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with 

G±(^t) = i. (5.5) 

According to the structure exhibited so far, we expect only even powers of A^ to 
occur in Gf n (A^). This is in accordance with /C and K? from ( f2.26| )-( p.27| ), as 
appearing in the Wheeler-DeWitt operator ( |2.15| ), being even in A and A^ . 

Due to our construction, the combinations (A^) n Gf n (A^) can be expected to be 
analytic at A^ = 0. One may however ignore the reasoning of Section 4 and treat 
( |5.4|) and (|5.5|) as an ansatz by its own (leaving the numbers and E n unspecified 
as well). This is the strategy we will pursue in what follows. The procecure is 
again to separate orders of a and to determine the solutions by some additional 
requirement. One might impose analyticity of (A*) n Gf n (A*) , but it turns out that 
a weaker condition does the job as well. We simply demand that 

(Hi) G^ n (A^) admits a Laurent series at A j = 0. 
In other words, Gf n (A^) may be expanded in integer (positive and negative) powers 
of A^ . Logically, this replaces the condition (i) as used in Section 4 (whereas an 
analogue of condition (ii) is no longer necessary). The technical point of condition 
(Hi) will be to exclude terms of the type ln(.4.t). 

The Wheeler-DeWitt equation in the Fock representation is of fourth order in 
derivatives with respect to A* . Although one might expect that this feature provides 
an additional complication, things actually become simpler. Writing a particular 
a) as JF(a, ^4+) |0) , we apply the Wheeler-DeWitt operator and thereafter divide 
the result by T(a, A^). This allows for a proper separation of orders of a. Proceeding 
iteratively, one encounters only first order differential equations of a very simple type. 
Also, the solutions Gf n (A^) turn out to be polynomials in positive and negative 
powers of (A*) 2 , hence are represented in terms of elementary functions. In other 
words, the Laurent series whose existence is required by condition (Hi) are actually 
finite. This is a great simplification as compared to the procedure of Section 4, 
where the coefficient functions F r ((f)) and G>(0) emerged as infinite series. 

The first non-trivial order 0(a) is of purely algebraic type and yields E n = (n + 
|)m, thus re-introducting the well-known eigenvalues of E. At 0(a°) we encounter 
the differential equation 
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the general solution of which is 
Gfn(^) = ± ^fe^ + Ki - ± (1 +p + 2q£)) ln(^) ± 3 {A \f , (5.7) 



where «i is an arbitrary integration constant. Condition (in) implies that the co- 
efficient of the ln(^) term must vanish. This is an equation for q^, the solution 
immediately turning out to be ( |5.2| ). Inserting all results obtained so far into the 
equation at 0(a _1 ), we obtain a differential equation of the type 

A^g'(A^ = p(A^, (5.8) 

where g(A^) stands for G2 n (A^) and p(A^) for an expression that has already been 
determined (up to the constant /%). Moreover, p(A^) contains only integer powers, 
ranging from (A^)~ 4 to {A^) A . The solution is thus a function of equal type, including 
an additive integration constant K2, and a ln(.4J) term whose coefficient turns out 
to be 

(-9E* - Am 2 El =f Am 2 E n ± S^m 2 ) . (5.9) 

Again, this term has to vanish. Thus K\ is fixed, and the complete solution for 
Gf n (A ] ) ^ 

GUA) = ± + IT* H ± 9 ± 4m ^) ± ■ ( 5 - 1Q ) 

8m{A>) z 8m/ v ' 8m 

This pattern persists at all orders. At 0(a~ 2 ) the differential equation for Gf n (A^) 
is again of the type Q5.8Q , with p(A^) being a finite sum of even integer orders of 
A\ and involving the constant k,2- The solution for Gf n (A^) thus consists of a finite 
sum of even integer orders of A\ an additive integration constant k 3 and a In (.At) 
term whose coefficient has to vanish, thus fixing K2, and so forth. We just display 
the highest and lowest order of the second coefficient function 

± f 9n(n-l)(n-2)(n-3) 9 t 4 

128m2(^t)4 + --- + 128m2^ j (5 ' Uj 

and note that in general Gf n (A^) contains contributions from (A^)~ 2r to (A^) 2r . 
Moreover, the structure is such that the combination Gf n (A^) = (A^) n G^ n {A^) is a 
polynomial. This is reflected by the coefficients n(n — 1) and n(n — l)(n — 2)(n — 3) 
in ( |5.10| )-( |5.11| ). Thus, despite the expansion in powers of a -1 , the regular nature of 
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the functional dependence on A* remains intact. As already mentioned above, such 
a feature is not at all generic for functions that are regular in two variables, but it 
serves here as part of the property singling out the wave functions fi„ . 

Since the oscillatory domain is bounded in a (and, moreover, exists only if n ^> 
m~ 2 ), there is no analogous expansion there. However, due to the WKB matching 
procedure as applied in Section 3, we define another set of exact wave functions 
by (|3.12j) (or, inversely by ( |3.13|) -( |3.14j )), where 9 n is still given by (|3.7p . In order 



to achieve the (approximate) oscillatory behaviour as in (|3.15| ), we have to define 



the normalization constants in ( |5.1| ) and ( |5.4| ) so as to give fl^ the asymptotic form 
(|310D - (EmD , hence 



This completely determines our wave functions fl^ and . Due to their structure 
(the appearance of the oscillator basis \n) at leading order) it is clear that a rather 
large set of solutions of the Wheeler-DeWitt equation may be expanded into them. 
By assuming 

oo oo 

ij> = 5>+ fi+ + k~ n~) = £(c+ + a' ~-) (5.13) 



n=0 n=0 



for some solution ip of the Wheeler-DeWitt equation, the k „ (or c„) may be com- 
puted by fixing some (large enough) a and expanding the initial data ip and d a ijj in 
terms of the at the fixed value of a (which should be possible if the bevaviour in 
is not too catastrophic). We will find a more convenient method later on, and the 
present argument just serves for the count of degrees of freedom contained in fl^ ■ 
Thus we treat {fi^} (or likewise {5^}) as a basis of the space M of wave functions, 
and a precise definition of which coefficients in ( |5.13| ) are allowed will be given later. 



The are, by construction, real (complex conjugation * being defined by its action 
in the position representation), while the transformation of the basis ( |3.12j ) implies 



/ ' — 'n 



(5.14) 



If n 3> m~ 2 and at the level of the WKB approximation at which fl3.15| ) is valid, this 
property is in accordance with and S~ representing an ensemble of collapsing 
and expanding universes, respectively. 

If one accepts the wave functions to play a distunguished role, one ends up 
with a distinguished decomposition of the space M of wave functions into the span 



26 



H + of and the span M~ of {H™}, hence M = M + © M~ . Wave functions 

are thus uniquely decomposed as ip = ip + + ip~ ■ It is important here to note that 
the differential structure of minisuperspace is not sufficient for identifying collapsing 
and expanding modes exactly. This is a major difference to the flat Klein Gordon 
equation (where a differential background structure, namely a timelike Killing vector 
field, enables one to define negative/positive frequency modes exactly in a Lorentz- 
invariant way) and it provides one of the most severe problems in constructing a 
consistent setup for quantum cosmology. However, if the analyticity structure of 
wave functions is accepted as a guiding principle in our model, we seem to have 
uniquely defined a decomposition of the space of wave functions which — in the 
WKB approximation — - is identified with collapsing and expanding (incoming and 
outgoing) modes. 

There is a formal relation between f2+ and Q~ that might provide a hint towards 
a deeper significance of the analytic structure we have exhibited so far. The Wheeler- 
DeWitt equation remains invariant under the substitution 

a — > ia m — > im , (5.15) 

leaving <p unchanged. Thus the quantities ma 3 and m~ 1 a (and by definition £, 
which involves m 1//2 a 3//2 ) remain unchanged as well, but we have a 2 — > —a 2 , and 
m 2 — > —m 2 . Under this substitution the wave functions fi+ and Q~, when written 
down in the large-a expansion, are almost perfectly transformed into each other. 
This includes the structure of the exponential prefactors as well as the normalization 
( |5~T2|) (cf. (151(5D-(PID and note that E 2 -> —E 2 while a/E n -> a/E n ), except for 



the prefactor a _1 / 2_p / 2 and the numerical factor ~ in fi~. (Also Hawking and Page 
have noted that this transformation carries their \l/^ p (a, 0) — which is our fig (a, <fi) 
- into the exponentially growing part of the no-boundary wave function as given by 
Page - - which is just our Qq (a, 0)). A similar relation exists between the Airy 
functions Ai(x) and Bi(x) in the expansion for large x, if x = a 4//3 is set. One could 
thus try to define M ± as the "eigenspaces" under a suitable substitution operation. 
However, this forces us to treat m as a variable rather than a fixed constant. We 
leave it open whether one would gain any new insight by doing so. 

The crucial question in exploiting the emergence of the spaces M ± for further 
developments of the subject of quantum cosmology is certainly whether the struc- 
ture showing up here carries over to the full (non-minisuperspace) Wheeler-DeWitt 
equation. 
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6 Scalar product and Hilbert spaces 

There is a natural (Klein Gordon type) scalar product associated with the Wheeler- 
DeWitt equation flU fll^l . If ip\ and ip2 are two solutions of the latter that are well- 
behaved for large <j), the expression 

poo / «-> \ 



(with f d a g = f d a g — (d a f) g) is independent of a on account of (|2.9|) . It defines 
an indefinite scalar product, the integrand being the a-component of a conserved 
current [pp . Due to its indefiniteness, it does not enable us to define a Hilbert space 
directly. 

When wave functions are expressed in the energy representation as introduced 
in Section 2, we find 

i 00 / <-> 

Wl, fa) = - - a P J2[ /n(°) d a 9n(a) + 
1 n=0 V 

Y a \/n(n + l) (f*^(a)g n+1 (a) - /* +1 (a)#„_i(a)) , (6.2) 

where ^i(a,£) = En/n(a)*n(0 ancl ^2(0, = En #n(a)*n(0- In tne Fock repre- 
sentation the scalar product reads 

% 3 
Q&i, ^2) = - - a p a| <9 a + - /C |^2, a) , (6.3) 

where the derivative 9„, when acting to the left, does not include the prefactor a v . 

Inserting the asymptotic expressions of our wave functions in the oscillatory 
domain, we find 

Q(p,f, nf) = Q(fi+ nj) = 5 rs (6.4) 

and 

Q(Ef, Zf) = ±6 rs Q(E+S;) = 0. (6.5) 

Using the asymptotic series for a —>■ 00, we find hints that these relations hold exactly 
for all r and s. For all those combinations in which the exponential a-dependent 
prefactor decreases, they are evident (Q being evaluated at arbitrarily large a). 
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For all other cases, we have used the first few terms of the series to check them. 
In the following we will assume that they hold. They also explain the particular 
normalization we have choosen (see ( |3.15| ) and ( |5.12p ). 

Let us consider now the decomposition M = H + @1H~ of the space of solutions of 
the Wheeler-DeWitt equation as induced by the wave functions H*. If any solution 



ip is expanded as in ( 5.13 ), the above normalization of the basis yields 

c± = ±Q(H± V). (6.6) 

Encoding the information contained in a wave function in terms of the numbers c^, 
we can talk about an "energy representation" in a sense quite more sophisticated 
than the notation of Section 2. Since the basis elements are solutions of the Wheeler- 
DeWitt equation, there is no dependence on additional variables, except for the 
"true" physical labels n and ±. 

Evidently, the scalar product Q is positive/negative definite on 1H + and M~ , 
respectively. Hence, admitting only wave functions if> whose coefficients satisfy 



oo 

£( 

n=0 



|c+| 2 + |c;| 2 )<oo (6.7) 



makes (H + ,Q) and (M~,—Q) two Hilbert spaces that may be used as a starting 
point for despriptions how to compute probabilities for observations. Admitting 
more general wave functions, one should still be able to compute relative probabili- 
ties. Given a solution if) as in ( |5. 13 ), the (relative) probabilities associated with the 



states are \Q(E^,ip)\ 2 . These numbers should be relevant when predictions for 
matter the energy contents are drawn. In the WKB-philosophy — which is relevant 
for quantitative predictions — this is quite clear for any basis being normalized 
as (|6.5|). (Any such basis gives rise to a decomposition into two Hilbert spaces). 
In Ref. 0] we have suggested a "minimal" interpretational scheme for quantum 
cosmology based on (IH, Q) as the only fundamental mathematically well-defined 
structure. However, here we suggest a different thing. If one accepts the basis wave 
functions as distinuished objects (not just as wave functions which describe a 
semiclassical ensemble of universes with approximately conserved matter energy), 
one is faced with a distinguished Hilbert space structure, based on asymptotic ana- 
lyticity properties. One might think of it to be introduced by some path-integral or 
to constitute a first principle by its own, and it constitutes an exact fundamental 
structure in addition to (IH,Q). If an analogous feature may be found in the full 
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(non-minisuperspace) well, it should be of some relevance to the fundamental 

conceptual problems of quantum cosmology. 



7 No-boundary wave function 

In order to show an example, we estimate the coefficients of the expansion (|5.13|) 



for the no-boundary wave function JOd ?/> nb of the Hawking model 0. It has been 
studied in great detail by Page ||. The classical no-boundary trajectories start from 
a point (do, 0o) on the zero potential line m 2 a 2 cj) 2 = 1 and are in the inflationary 
domain ma|0| 3> 1, |0| 3> 1 given by 

1 exp^(^-0 2 )Uexp^(^-0 2 )V (7.1) 



This expression is valid as long as 1 |0| <C \<f> — l/(30 o )|. Since there is only 
one no-boundary trajectory hitting any point (a, 0) in the inflationary domain, we 
associate with this point the initial value 0o of the latter, thus turning O into a 
function of a and 0. The action of this congruence of trajectories is 

s = - w ( m2<iV - »r + 5 * - • < ? ' 2 > 

the no-boundary wave function in the inflationary domain being (thereby general- 
izing Page's expression to arbitrary p) 

^ NB (a, 0) « a- p/2 - 1 A((f> ) cos(S(a, 0)) (7.3) 



with 



A(cf> ) = 




and 0o now interpreted as a function of a and 0. Its basic structure is that it is a 
product of the rapidly oscillating WKB-type function cos(S') with a slowly varying 
prefactor. (The prefactor A(0 O ) is in fact to some extent arbitrary. In the WKB- 
approximation the Wheeler-DeWitt implies that it is constant along the classical 
trajectories. The expression ( |7.4|) corresponds to the solution arising from the no- 
boundary proposal for the Euclidean path-integral). 
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The no-boundary trajectory with initial value 0o leaves the inflationary domain 
(i.e. attains \<f>\ ~ 1) at a min rs m -1 ^! -2 ' 3 ^ ' 2 , and subsequently undergoes the 
matter dominated era in which oscillates and the matter energy E is approximately 
conserved. Since, on the other hand, a min rs (n/m) 1 / 3 , we find 0q ~ | In (nm 2 ), which 
makes 0o in (|7.3| ) effectively a function of n. Since |0 O | 3> 1, we have n 3> m~ 2 , and 
thus a non-trivial classical domain. The prefactor (|7.4| ) then becomes 

as far as the contribution of trajectories representing universes with matter energy 
E n is concerned. 

We will choose an indirect way to estimate the magnitude of the coefficients 
when ?/> NB is expanded into 5^ as in (|5.13| ). Since we do not know precisely the be- 
haviour of V^nb^, 0) in just those regions in which we know S*(a, 0), the information 




necessary to perform the integration (|6.6|) is not easily accessible. We may instead 
first evaluate the relative probability distribution for trajectories labelled by O m 
the position representation. Since Q(ip,ip) = for real if), we consider the incom- 
ing/outgoing projections i/j^ b ~ |a _p / 2_1 Ae TlS . The expressions Q("0nb; ^nb)> wnen 
computed according to ( |6.1|) , turn out to be integrals over the measure ±|mA 2 \<f>\d<j) 
(the additional a-dependence cancelling, as it should). Using \(f>\d(f> ~ \(f> \d(f>o (at 
constant a), this becomes ±P(0o) d<f>o (± the relative probability for finding the 
universe represented by a trajectory in the interval between 0o and 0o + <i0o), where 

777 

P(0o) = j |0o| ^(0o) • (7.6) 

This is the standard procedure of evaluating probabilities for WKB-type wave func- 
tions based on the conserved current. (An alternative interpretation || predicts 
probabilities that differ from these by the amount of proper time spent by trajec- 
tories in a domain of minisuperspace). By including an additional factor 2, we 
may restrict O to be positive. Since in this case O an d n are related uniquely 
by O « \\J2 ln(nm 2 ), we may compute d^jdn « (9ri0 o ) _1 . Due to the symme- 
try ip NB (a, — 0) = -0 NB (a, 0) we know that = for odd n. Supposing that |c^| 
for even n may be approximated by continuous functions, we set dn = 2 and find 
that the relative probability for the universe to have energy quantum number n is 
given by P n ~ 4 (9n0o) _1 P(0o)- On the other hand, the expressions Q(4>nbi V'nb) 
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are given by ± J2 n | | =l= , and the relative probability to find the universe containing 
matter energy E n in the contracting/expanding mode is = ~P„ = |c^| 2 for even 
n (and zero for odd n). We thus identify for even n and n 3> m~ 2 



, ,„ 1 2P(6 ) 1 [m A 

^ n ; 3 V ™</>o 3 V 2n 

1 1 



(y / 6 + 2(exp( — (7.7) 

^21n(nm 2 ) V V Pl 2m 2 ln(nm 2 ) ; )) n ^ > 

where \K n \ = 1. The first equality is exact, it stems from the fact that is 
real. (Note that due to the reality of V'nb, both ± modes are of equal probability, 
although this number is infinite). 

The density P((f>o) as well as the probabilities P n display the well-known problems 
for the no-boundary wave function to predict sufficient inflation (i.e. a sufficienty 
large universe; see Refs. jpT|[|18|). Due to the smallness of m, the exponentials 
prefer small <fio and n, and it is only the flat behaviour of P{4>q) for <ft — > oo (or the 
dominant behaviour P n ~ rT 1 asn^ oo) that allow for large universes. However, 
trajectories with m|0o|~l correspond to classical universes above Planckian densities 
after nucleation, and it is not clear whether they should contribute (cf. Ref. |J). 
These trajectories correspond to n>,m~ 2 exp(^2-) ~ exp(5 x 10 12 ) (and thus a max 
being given roughly by the same number, or, expressed in "true" units as displayed 
in the Introduction, a max rs exp(5 x 10 12 ) centimeters, or light years, or present 
Hubble scales, which makes no big difference due to the huge value of this number). 



As a consequence of (|7.7|), we find in the energy representation |t/>nb,o) = 
J2nfn{a)\n) that the oscillator components (for even n 3> m~ 2 and in the range 



(n/m) 



1/3 < 



0«Oj 



2mn) are given by 



fn(a) 



2 a -p/2-l/4 t /™ 



.4, 



2n $2EZ 



cos ( = {2E n ) 1/2 a 3/2 + 5 n 



(7.8) 



Here, we have set K n = e tSn . Due to the large value of n one may of course set 
E n w mn. The components for odd n vanish. For a>,a max the / n (a) contain the 
exponential e ±a2 / 2 ^ E " a terms of (]3.10|) -( 3.11|) . An analogous large-a behaviour is 
expected to apply for the component functions with n<,m~ 2 (which do not contribute 
to classical universes) as well. 
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The expression ( |7.8| ) may heutistically be checked by using the real part of 
(|2.32|) with k = |m 1 / 2 a 3//2 as an approximation for cos(S'). Invoking \l/ n (x) ~ 
\E f n (0) cos(x\/2n) for small x and even n, and \l/n(0) ~ (— ) n / 2 (2/n) 1 / 4 7r -1 / 2 for large 
even n (which follows from Stirling's formula for n\), a behaviour roughly similar 
to (|7.8| ) is recovered but without the phases 5 n and an additional factor | in the 
Cosine. This reflects our lack of knowledge about the details of V'nb^, 4>) m the do- 
main where S^(a, 0) is known, and vice versa. It would be interesting to study this 
problem in more detail, and to find an estimate for the 5 n . In case these numbers 
vary rapidly with n, the Cosine in ( [7.8| ) would introduce a chaotic type behaviour 
of the oscillator excitations. 

In their work on wormholes Hawking and Page have assumed that the no- 
boundary wave function (which increases as e a I 2 ) can be expanded in terms of their 
wave functions \l/^ p (a, 0), which are our fi~(a,0) and decrease as e~ a / 2 . In our 
language, they have assumed to form a basis, in which case one would expect 

the expansion coefficients showing a tremendous increase with n. In contrast, in 
our framework, these states are only half of a basis. Note that linear independence 
relies on the precise definition of a vector space. (For example, Sine and Cosine 
functions can or cannot be expanded into each other, depending on the interval 
on which they are considered). However, by just counting degrees of freedom, the 
existence of a framework in which ?/>nb may be expanded in terms of is very 

unlikely. (This would in fact imply that f2+ may be expanded in terms of Q~, and 
the structure provided by the scalar product Q would break down). The explicit 
computation of the expansion coefficients k„ with respect to as defined by ( |5.13|) 
seems to be difficult at the present status of our knowledge. For c+ = (c~)*, we 
find k+ + ik~ = 2c+e 40n with O n from (|3.7|). Our estimate (|7.7|), together with 
K n = e tSn , leads to a real expression multiplied by e^ e " +s "\ Since we do not know 
5 n , we cannot compute the real and imaginary part of this phase factor directly. On 
the other hand, an estimate for 6 n might be accessible by using further information 
about ipNB- The mere fact that it containts an increasing e a 7/2 contribution implies 
that at least some O n + 5 n are different from | (modulo 2n). 

There is another interesting question related with the work of Hawking and Page. 
They assume that the regular superpositions of {^~} — due to their exponentially 
dampted nature — provide quantum wormhole states. On the other hand, the sec- 
ond half of the basis appears on quite an equal footing here: both types of 
states describe universes with matter energy E n . The only possible essential dif- 
ference concerns observations when the universe is near its classical turning point 



33 



o- ~ flmax • It is in particular the states Q~ that do not seem to provide problems 
there. How does this fact relate to the interpretation of certain superpositions of 
these states as wormholes? Usually, the exponentially damped behaviour e~ a I 2 is 
associated with wormhole states by definition. However, we do not have a satisfac- 
tory interpretation of the exponentially increasing behaviour e a I 2 (nor is it required 
by definition for any wave function, but just accepted as a grain of salt rather than 
a desired property when it emerges). This is a certain asymmetry (at least as long 
as no path-integral arguments are invoked) that might point towards a theoretical 
lack in our understanding of quantum cosmology. We are not able to give an answer 
to this question, but it seems worth pursuing it. 

The predictive power contained in the coefficients (|7.7f ) is — as far as practical 
quantitative features are concerned — equal to the results of the common WKB- 
philosophy. If, however, the decomposition M = 1H + © M (with IH being defined 
by the normalization condition 0fc>.7| )) is regarded as a distinguished one, we may 
embed V'nb into a mathematically well-defined underlying structure. We have two 
fundamental Hilbert spaces (M ± ,±Q) and a wave function for which Q(H,-?/;nb) is 
finite for any element H G M, and it is just these numbers in which all physical 
information about ?/>nb is encoded. The wave function i/j^b is not an element of M. 
Since both projections i/j^ b onto H ± have infinite norm Q(V'nB) ^nb)i they are of 
distributional character (similar to momentum eigenstates in conventional quantum 
mechanics) . 



8 Concluding remarks 

The behaviour of the real wave functions f^(a, 0) for large a provides a relation 
between the expansion ( |5.13| ) of some ip into these (the coefficients k„ being trivially 
connected to by ( |3.12| )-( |3.14j )) and the question of boundedness of ip- Due to 
the factors e a / 2 , all f2+(a, 0) are unbounded. Hence, a given wave function ift(a,(j)) 
seems to be bounded away from a = (i.e. \ip(a,(fi)\ < K < oo in any domain 
a > ax) if and only if = for all n. (Note that this is in contradiction with 
the remark in Ref. concerning the expectation that V'NBj which is not bounded, 
may be expanded in terms of \E^ P = fi^). It is not entirely clear to what extent 
the unboundedness of a wave function causes interpretational problems (e.g. near 
the classical turning point). Some authors consider boundedness as a condition 



necessary for interpretation [19|]B and talk about a "final condition" for the wave 
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function. Such approaches could provide an additional justification for expecting 
the limit a — > oo to play a conceptually fundamental role. At the technical level we 
do not know whether the degrees of freedom contained in the coefficients k~ may 
be arranged so as to cancel the expected singular behaviour of Q~ (a, 0) for small 
a and make up a strictly bounded solution (although Hawking and Page, when 
constructing approximate wormhole states, provide a hint that this is possible). In 
the case of the tunnelling wave function emerging from the outgoing mode proposal 
20flpl], the boundedness of ?/>t is usually considered part of the definition. If such 



functions exist at all, this would immediately imply that ipT is a superposition of the 
Q~ alone. It might be worth thinking about whether a possible relation between ip^ 
and the wormhole context in which Hawking and Page considered the wave functions 
^ p = Q~ sheds some new light on the conceptual questions. 

Let us close this article by adding some general speculations. If the structure 
encountered in the Hawking model carries over to some more sophisticated (prefer- 
ably non-minisuperspace) model one would apply WKB-techniques in combination 
with decoherence arguments (traces in the Hilbert spaces M ) in order to identify 
states with physical observables and to recover the standard laws of physics. The 
Q-product should boil down to plus or minus the standard scalar product of quan- 
tum mechanics fl | [ P^j . The existence of two separate Hilbert spaces may be a hint 



that the ± sectors decouple from any observational point of view: given a wave 
function, one is either in the + or in the — sector, no experience of a superposition 
is possible. One could call this a super-selection rule. Only within these sectors 
Hilbert space techniques apply, and the actual non-experience of various other su- 
perpositions is delegated to decoherence. The ultimate object to describe experience 
would be a reduced density matrix, as evaluated by standard Hilbert space methods, 
hence within a completely well-defined framework. It remains to be seen whether 
such an interpretation is still possible when observations near the turning point are 
concerned. There, one would heuristically expect to undergo a "transition" from a 
reduced density matrix belonging to the — sector to one belonging to the + sector. 

Since the situation is a bit reminiscent of the one-particle Hilbert spaces of neg- 
ative/positive frequency modes as emerging from the flat Klein Gordon equation, 
a further possible direction to pursue is to envisage a third-quantization |2^j in 
terms of the preferred decomposition (cf. Ref. fl2fl). 
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